Broken symmetries and directed collective energy transport 
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We study the appearance of directed energy current in homogeneous spatially extended systems 
coupled to a heat bath in the presence of an external ac field E{t). The systems are described by 
nonlinear field equations. By making use of a symmetry analysis we predict the right choice of 
E{t) and obtain directed energy transport for systems with a nonzero topological charge Q. We 
demonstrate that the symmetry properties of motion of topological solitons (kinks and antikinks) 
are equivalent to the ones for the energy current. Numerical simulations confirm the predictions of 
the symmetry analysis and, moreover, show that the directed energy current drastically increases 
as the dissipation parameter a reduces. Our results generalize recent rigorous theories of currents 
generated by broken time-space symmetries to the case of interacting many-particle systems. 
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The idea of rectifying energy transport with the help of 
fluctuations has been discussed for several years in con- 
nection with molecular motors and other nonequilibrium 
properties of biological systems ^ , electrical currents in 
superlattices voltages in Josephson junction coupled 
systems to name a few. The fluctuations have zero 
mean value, i. e. the dc component is absent. The reduc- 
tion of the underlying problem, namely directed energy 
transport, to a particle moving in a space-periodic but 
asymmetric (ratchet) potential allowed to study the re- 
sulting directed current in great details [Q (for a recent 
review 

see Q). A recently elaborated symmetry approach to 
this problem established a clear relationship between di- 
rected currents and broken space-time symmetries [^0- 
The essential step was to separate the unavoidable cor- 
relations in the fluctuations from the uncorrelated ones. 
This is easily obtained by replacing the fluctuations as a 
superposition of ac driving fields and uncorrelated white 
noise. An important consequence is that the symmetries 
may be broken either by violating the reflection sym- 
metry of the potential in space or by violating the shift 
symmetry of the ac fields. Thus, a particle may display 
a directed motion also in the case of a space-symmetric 
potential. Another interesting result is the persistence 
of directed currents in the Hamiltonian limit of systems 
exposed to ac fields but decoupled from the heat bath 



A very important question is whether the symmetry 
approach can be generalized to the case of interacting 
many-particle systems (for related studies see P,p|-p^). 
In this Letter we focus on nonlinear partial differential 
equations, which can be considered as classical analogues 
of quantum models of interacting particles |11|. We show 
that such systems allow for a directed energy transport 
when being driven by a proper combination of ac forces. 
The success of such a generalization will not only un- 



derline the general validity of the symmetry analysis, 
but will also allow for a systematic symmetry analysis 
of nonadiabatic nonlinear response functionals for many- 
body theories. 

We study the properties of nonlinear Klein-Gordon 
equations for a scalar field tp which depends both on a 
spatial coordinate x and time t: 
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where cq is the limiting propagation speed of small am- 
plitude plane waves, a is the dissipation parameter de- 
termining the inverse relaxation time in the system, and 
= df /dz. The ac field E{t) has zero mean and pe- 
riod T. The Gaussian white noise ^ is characterized 
by the standard correlation function {^{t,x)^{t' ,x')) = 
2^6{x — x')d{t — t') where (3 is the inverse temperature. 
The potential U{z) is assumed to have several identi- 
cal minima for z = Zi, i = 1,2,.... We choose fixed 
boundary conditions: ip{x — oo) = zi, ip{x +oo) = 
z„i. These boundary conditions determine the topolog- 
ical charge Q = m — I which will be of importance in 
the following. By making use of the standard continu- 
ity equation p^t + j,x = , where the energy density 
p — {if \ -\- c^tp\)/2 + U{(p), we obtain the following ex- 
pression for the energy current J: 



J{t) 



rdx . 



(2) 



Next we turn to the symmetry analysis of (|ljy). It is 
carried out in the absence of ^ in (|l|) , as it does not affect 
the symmetries. We note that J changes its sign under 
the transformation x — > —x. Eq. (|l|) is invariant under 
the space inversion. However the only topological charge 
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invariant under this transformation is Q = 0. Conse- 
quently we conclude that for zero topological charge, 
(3 = 0, the time-averaged energy current (J) vanishes 
exactly. 

In order to study the possibility of the appearance of 
a nonzero energy current in the presence of a nonzero 
topological charge, Q 7^ 0, we note that the topological 
charge Q 7^ is invariant under the combined symmetry 
operation x ^ —x and ip —ip + zi+Zm- The same 
symmetry operation changes the sign of the energy cur- 
rent, and all we need is to ensure that the equation (|^) 
remains invariant. This is now only possible if the ac 
field possesses shift symmetry and the potential U{z) is 
symmetric around z = {zi + Zm)/'2. Thus, we arrive at 
the symmetry operation 
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that leaves the equation ([^) invariant and changes the 
sign of the energy current as the constraints 



E{t) = -E{t + T/2) 
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hold. In such a case we may conclude that the directed 
energy current vanishes Violating (jj) we loose the 
symmetry and consequently may expect a nonzero 
energy current provided no further hidden symmetries 
are overlooked. 

Most interesting is that the conditions can be vio- 
lated by choosing a function E(t) which is not shift sym- 
metric. A simple choice is 



E{t) = El cos ujt + E2 cos{2Lut + A) , 



(5) 



with El ^ and E2 7^ 0. In particular, such an ac field 
allows to generate a nonzero energy current in the well 
known sine-Gordon equation with 



U{z) 



(6) 



which we will consider in the following. At the same time 
we stress here that any other choice of the ac drive E{t) 
which is not shift symmetric (e.g. a corresponding se- 
quence of pulses) will do the job, as well as any other 
potential U{z) which meets the above described condi- 
tions. 

Of special interest is that in the Hamiltonian limit as 
a ~> 0, time reversal symmetry is recovered: the symme- 
try operation t —t leaves invariant, provided E{t) 
is symmetric: E{t) = E{—t). As time reversal always 
changes the sign of the energy current (||), zero average 
energy current will be a consequence. Thus, for the un- 
derdamped case a <C 1 we expect a change of sign of the 
energy current upon varying A, because for the partic- 
ular values of A = 0, tt time reversal is approximately 
restored and the directed energy current disappears. 



As discussed above, a necessary condition for a directed 
energy current is a nonzero topological charge Q. It im- 
plies the presence of an excess of kinks over antikinks or 
vice versa. The simplest case is Q = 1 as a single kink 
is present on average in the system. As a kink possesses 
a nonzero rest energy Ek = 8co, it is natural to expect 
that the nonzero energy current is generated by a directed 
motion of such a kink. Indeed, assuming that the kink 
position X and velocity V are described by expressions 



xip^xtdx . 

(7) 



we observe that the symmetry operation changes the 
sign of V. Thus, we find that if the energy current van- 
ishes by symmetry, the same is true for the average veloc- 
ity of a kink. This shows the intimate connection between 
a directed energy current and coherent kink excitations as 
carriers of such an energy. In [l0| a directed kink motion 
was obtained numerically for the first time for the case 
(^,(||) together with a study of its complex dynamics. 
Note that in a directed kink transport has been 

obtained by breaking the reflection symmetry of U{z) for 
^2 = 0. 

If the temperature and thus the amplitude of ^ is large 
enough, additional kink-antikink pairs are excited in the 
course of evolution. Since the space inversion (x — > —x) 
does not change eq. (|l|), transforms Q ^ —Q and kinks 
into antikinks with inverted velocities, the average ve- 
locity of a kink-antikink pair is zero. Consequently the 
generation of kink-antikink pairs does not affect the net 
energy current of the system predefined by the value of 
topological charge Q. 

With these analytical results we turn to a numerical 
investigation. We used a standard discretization scheme 
and integrated the following coupled differential equa- 
tions: 

4>n - cl (0„-l - 2(j)n + (jjn+l) + "^n + ^T" = /(*) • (8) 

dcpn 

Here n is an integer representing the discretization in 
space. We choose the following parameters: /3 = 100, 
4 = 10, El = 0.2, E2 = 0.2, uj = 0.1 and Q = 1. This 
particular value of cq ensures that the effects of discrete- 
ness on the kink motion may be neglected EbI. 
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FIG. 1. Dependence of the time-averaged energy current 
(J) on the phase shift A. Solid line - a = 0.05, dotted line - 
a = 0.01, dashed-dotted line - a = 0.2. Note that the latter 
case a = 0.2 is scaled by a factor of 5, so the real values of 
(J) are five times less than they appear in the plot. 

The time averaged energy current (J) as a function of 
the phase shift A for three different values of dissipation 
parameter a = 0.01,0.05,0.2 is shown in Fig. 1. We ob- 
tain a strong dependence of the directed energy current 
(J) on A with sign changes, as expected. Furthermore 
the maximum values of the energy current increase sub- 
stantially as the dissipation parameter a reduces, similar 
to the case of a single particle moving in a space-periodic 
potential 0. Especially we find an increase by a factor 
of 50 when the damping a decreases from 0.2 to 0.01. 
And finally we clearly observe that as the parameter a 
reduces, the zeroes of the curves tend to the positions 
determined by time reversal symmetry, A = 0, tt. We 
computed in a similar way the time averages of Y as 
defined in (^, and obtained identical results upon vari- 
ation of A and a. It shows that indeed kink motion is 
responsible for the energy transfer. 

In these simulations the amplitude of the white noise 
^ was small enough to prevent an excitation of a kink- 
antikink pair during the simulation time. Nevertheless a 
very long simulation should also show up with a spon- 
taneous creation of such pairs. In order to see that this 
event does not influence the symmetry analysis and es- 
pecially does not change the obtained values of the en- 
ergy current (cf. above considerations) we proceed with 
a comparison. In Fig. 2 we plot the dependence of the 
time-averaged energy current (J) on the time of averag- 
ing for a — 0.01 and A = 7r/2, as a single kink is present 
in the system. In the same figure we plot the result for 
the case of additionally exciting a kink-antikink pair Jl7[ . 
We indeed observe that the value for the time-averaged 
energy current does not depend on the number of ad- 
ditional kink-antikink pairs. As the simulation time in- 
creases even further, soliton collisions eventually lead to 
the dissipation of the pair, with one single kink remaining 
in the system. 




FIG. 2. Dependence of the time-averaged energy current 
(J) on the averaging time t for A = 7r/2 and a = 0.01 for the 
two cases: a single kink in the system (solid line), a single kink 
and an additional kink-antikink pair in the system (dotted 
line). The insets (a) and (b) show snapshots of the field ip as 
a function of x at some time during the simulation for these 
two cases. Arrows indicate the direction of time-averaged 
motion of the kink(s) and the antikink. 

Let us discuss possible mechanisms of energy current 
rectification in spatially homogeneous extended systems 
through directed kink motion. First we notice that in 
many cases the kink motion may be mapped to the clas- 
sical mechanics of a macroscopic particle iQ. However, 
a free particle that is subject to a spatially homogeneous 
ac force does not display a directed motion at all. One 
possible source of the observed directed current can be 
relativistic effects which are inherent to the considered 
nonlinear Klein-Gordon field equations. They can be 
explicitly obtained in some cases using projection tech- 
niques to derive effective equations of motion for the kink 
In such a regime the kink motion is mapped 
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onto the one of a free relativistic particle in the presence 
of damping and the ac field E{t). It is straightforward to 
see that this case leads to zero rectification for the limit of 
small damping, in sharp contrast with our findings. Thus 
we argue that the origin of the observed strong rectifica- 
tion in the underdamped limit is due to the nonadiabatic 
excitation of internal kink modes and their interaction 
with the translational kink motion 



|io|,yii9|. 

It is worthwhile to mention a number of other ways 
to generate a nonzero energy current. One case is to 
consider strongly discrete systems where the appearance 
of a nonnegligible Peierls-Nabarro potential can gener- 
ate a net kink motion through the lattice. Also all sorts 
of perturbations in the form of spatial modulations may 
act similar to discretization effects pO| , p]] |. Yet another 
way of breaking the symmetry has been analyzed in p^ ] 
where in addition to an ac force E{t) = cost the poten- 
tial U{z,t) — {z^ — 1)^ -I- sin t was assumed to be time- 
dependent. The application of our symmetry analysis 
yields that the symmetries (^ are violated and directed 
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energy transport occurs. It is also interesting to recall 
that a number of publications have been dealing with di- 
rected kink motion in asymmetric double well potentials 
| p5| , p3t . The asymmetry was chosen in a way to entropi- 
cally favour one equilibrium over the other one (still keep- 
ing the energy values of the equilibria identical). Even 
Gaussian white noise leads then to the directed motion 
of a kink. However additional ac fields have been shown 
to counterbalance the entropic force. 

While it is certain that our results may be applied to 
many different physical situations where the sine-Gordon 
equation is replaced by other Klein-Gordon equations 
with topological excitations, we stress that a particu- 
lar system, namely an annular long Josephson junction 
with a trapped fluxon, is perhaps the ideal case to ver- 
ify our results. This system is described precisely by (|l|) 
|^,||j25|. While in experimental results for E2 ^ 
have been presented, all one needs is to add the second 
harmonics E2 ^ and vary the phase shift A. The ana- 
logue of the ac field is here an ac homogeneous microwave 
radiation which traverses the junction. Under these con- 
ditions a trapped fluxon displays a directed motion, and 
correspondingly a dc voltage drop across the junction ap- 
pears. This effect is the analogue of the energy current 
discussed here. Notice here that the dc voltage drop os- 
cillates with the phase shift A and substantially increases 
in the "strong" underdamped Hmit, a ^ iv. 

We conclude this work with stating that the use of 
symmetries allowed to obtain a relationship between a 
directed energy current and directed kink motion. This 
may be of further interest as it shows an intimate con- 
nection between symmetry breakings and the relevance of 
coherent excitations and their properties in many body 
theories. Finally we note that nonlinear Klein-Gordon 
systems which do not admit kink or antikink excitations 
(zero topological charge) do not allow for directed energy 
currents induced by ac forces, and that the above consid- 
erations can be used without change for spatially discrete 
systems. 
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